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Abstract. The BL Lacertae object AO 0235+16 is well known for its extreme optical and radio variability. New 
optical and radio data have been collected in the last four years by a wide international collaboration, which 
confirm the intense activity of this source: on the long term, overall variations of 5 mag in the R band and up to a 
factor 18 in the radio fluxes were detected, while short-term variability up to 0.5 mag in a few hours and 1.3 mag 
in one day was observed in the optical band. The optical data also include the results of the Whole Earth Blazar 
Telescope (WEBT) first-light campaign organized in November 1997, involving a dozen optical observatories. The 
optical spectrum is observed to basically steepen when the source gets fainter. We have investigated the existence 
of typical variability time scales and of possible correlations between the optical and radio emissions by means of 
visual inspection and Discrete Correlation Function (DCF) analysis. On the long term, the autocorrelation function 
of the optical data shows a double-peaked maximum at 4100-4200 days (11.2-11.5 years), while a double-peaked 
maximum at 3900-4200 days (10.7-11.5 years) is visible in the radio autocorrelation functions. The existence 
of this similar characteristic time scale of variability in the two bands is by itself an indication of optical-radio 
correlation. A further analysis by means of Discrete Fourier Transform (DFT) technique and folded light curves 
reveals that the major radio outbursts repeat quasi-regularly with a periodicity of ~ 5.7 years, i.e. half the above 
time scale. This period is also in agreement with the occurrence of some of the major optical outbursts, but not 
all of them. Visual inspection and DCF analysis of the optical and radio light curves then reveal that in some 
cases optical outbursts seem to be simultaneous with radio ones, but in other cases they lead the radio events. 
Moreover, a deep inspection of the radio light curves suggests that in at least two occasions (the 1992-1993 and 
1998 outbursts) flux variations at the higher frequencies may have led those at the lower ones. 
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1. Introduction 

BL Lacertae objects, together with quasars with flat ra- 
dio spectrum, belong to the class of active galactic nuclei 
called blazars. These are objects characterized by extreme 
emission variability at all wavelengths, noticeable emission 
at the 7-ray energies, important polarization, superlumi- 
nal motion, brightness temperatures by far exceeding the 
Compton limit (see e.g. Urry 1999 for a review). 

The general scenario for the blazar emissivity has long 
been set: it foresees a plasma jet coming out from a su- 
permassive, rotating black hole surrounded by an accre- 
tion disk. In the jet, which is closely aligned with our 
line of sight, relativistic electrons produce soft photons 
(from radio to the UV-X-rays) through synchrotron emis- 
sion, and high-energy photons (up to the TeV energies) via 
inverse Compton scattering. However, theoretical models 
differ in the jet structure and physics, and in the origin 
of the seed photons that are energized to the 7-ray en- 
ergies. Constraints on the theoretical models can come 
from dense optical and radio monitoring, and simultane- 
ous observations throughout all the electromagnetic spec- 
trum, from the radio to the 7-ray band. Indeed, multi- 



wavelength variability studies (see Ulrich et al. 1997 for 
a review) can give information on the compactness of the 
emitting regions, and verify the existence of correlations 
and time delays among the emissions in different bands. 
This in turn can shed light on the location of the emitting 
regions in the jet and on the nature of the seed photons 
that are comptonized. International collaborations among 
optical astronomers such as the OJ-94 Project and the 
Whole Earth Blazar Telescope (WEBT) were born in the 
last years to make the observational effort more efficient. 

In this paper we present the results of a wide collab- 
oration aimed to study the optical and radio variability 
of the BL Lac object AO 0235+16. Partial and ve ry pre - 
liminary results were presented in Villata et al. ( 1999|) , 
where the light curve in the R band from January 1993 
to January 1999 was shown and radio-optical correlations 
during that period analyzed by means of the Discrete 
Correlation Function (DCF). 

The source AO 0235+16 is a well known blazar at z — 
0.94 that exhibits spectacular emission variations on both 
short (~ day) and long (months, years) time scales. 

It was classified as a BL Lac object by Spinrad & Smith 
(1975), who observed a ~ 2 mag variation. No ticeab le op- 
tical o utbur sts were descri bed by Rieke et al. ( 1976), Pic a 
et al. (|198C| ), Webb et al. (|1988| ), Webb & Smith (|l989p , 



Webb et al. ( 2000 ). From the data reported in the liter- 
ature, an overall brightness variation of more than 5 mag 
can be seen, which makes AO 0235+16 one of the most 
interesting sources for variability studies. 

As for the short-term optical behaviour, intraday vari- 
ability was o bserved in a nu mber of occasions (Heid t & 
Wagner |l996| ; Noble & Miller |l996t Romero et al. |2000| ) . In 



particular, changes of 0.5 mag were measured by Romero 
et al. (200C) within a single night, and variations of more 
than 1 mag in about 24 hours. 

Many optical data were taken in the period 1993- 
1996 by the astronomers participating in the interna- 
tional collaboration called OJ-94 Project, who observed 
AO 0235+16 as a "complementary" object, in addition 
to OJ 287 and 3C 66A. The results of their monitoring 
were published in Takalo et al. (1998). Some of the groups 
belonging to the OJ-94 Project have then continued to 
observe AO 0235+16. 

A huge monitoring effort has been pursued in the radio 
band by the University of Michigan Radio Observatory 
(UMRAO), in the USA, and by the Metsahovi Radio 
Observatory, in Finland. Indeed, this object is of extreme 
interest also in the radio band. Long-term radio light 
curves at frequencies above 1 GHz have shown that it 
exhibits frequent, relatively well-resolved, high-amplitude 
events (O'Dell et al. |1988|; Clements et al. |l995[). This 
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makes AO 0235+16 particularly suitable for studies of cor- 
related activity. Moreover, several recent papers have ar- 
gued that unusually high Doppler factors/large Lorentz 
factors may be present: Jorstad et al. ( [2001 ) infer a 
Lorentz factor greater than 45 based on 43 GHz proper 
motions, using VLBA observations during 1995 and 1996; 
a Doppler factor for the bulk flow of order 100 was in- 
ferred from radio short-term variability in an earlier pa- 
per by Kraus et al. ( 1999[ ). Doppler factors larger than 
80 (and Lorentz factors T > 60) were found by Fujisawa 
et al. ( [1999 ) by applying the inverse-Compton effect and 
equipartition models to VLBI observations at 22 GHz. 
Recen t VSOP observations of AO 0235+16 by Frey et al. 
( |2000| ) have led to an estimate of T b > 5.8 x 10 13 K for the 
rest-frame brightness temperature of the core, which is the 
highest value measured with VSOP to date and implies a 
Doppler factor of ~ 100. 

Optical and radio data on AO 0235+16 taken in 
the last four years (1996-2000) by a wide international 
collaboration are presented in Sect. 2, and inserted in 
the source long-term light curves, dating back from the 
mid Seventies. The results of the Whole Earth Blazar 
Telescope (WEBT) first-light campaign on AO 0235+16 
in autumn 1997 are also included. Optical spectra are de- 
rived and spectral changes discussed in the same section. 
Radio light curves are shown in Sect. 3, where the be- 
haviour of the radio flux is analyzed side by side with 
the optical one. Statistical analysis is presented in Sect. 4: 
autocorrelation function and Discrete Fourier Transform 
(DFT) methods are applied to search for characteristic 
time scales of variability in both the radio and optical 
domains; DCF analysis is then performed to check the ex- 
istence of radio-optical correlations. A final discussion is 
contained in Sect. 5. 



2. Optical behaviour of AO 0235+16 

The long-term monitoring optical light curves of AO 
0235+16 in the Johnson's UBV and Cousins' RI bands 
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Observatory 


Nation 


Long, (deg) 


Lat. (deg) 


Tel. diam. (m) 


Filters 


Symbol" 


MSSSO 


Australia 


+ 149.066 


-31.277 


1.0 


R 


filled triangles 


Xinglong 


China 


+117.575 


+40.393 


0.6 


R 


filled squares 


Vainu Bappu 


India 


+78.83 


+12.57 


2.0 


R 


asterisks 


Abastumani 


Georgia 


+42.80 


+41.80 


0.7 


R 


open circles 


Tuorla 


Finland 


+22.17 


+60.27 


1.03 


V 




Roma 


Italy 


+12.70 


+41.80 


0.5, 0.7 


BVRI 


open triangles 


Perugia 


Italy 


+12.372 


+43.112 


0.4 


VRI 


open diamonds 


Torino 


Italy 


+7.775 


+45.038 


1.05 


BVRI 


open squares 


Calar Alto 


Spain 


-2.546 


+37.224 


2.2 


R 


filled circles 


Teide 


Spain 


-16.5 


+28.16 


0.82 


BVRI 


crosses 


NOT 


Spain 


-17.88 


+28.75 


2.56 


BVRI 


asterisks 


Foggy Bottom 


USA 


-74.25 


+43.24 


0.4 


R 


filled stars 


Climenhaga 


Canada 


-123.22 


+48.25 


0.5 


R 


filled diamonds 



a The symbols are used to distinguish the R band data plotted in Figs. 3-7 



are shown in Fig. 1. As generally happens, the best sam- 
pled band in the past is the B one, while the use of CCD 
cameras has led people to observe chiefly in the R band. 
The data plotted before the vertical line at JD = 2450300 
were derived from the literature (see references in the cap- 
tion); those after are new observations presented in this 
paper. Thirteen optical observatories all over the world 
were involved in this project; their names, nationality, 
position (longitude and latitude), diameter of the tele- 
scope and filters used are indicated in Table 1. Seven of 
them (MSSSO, Xinglong, Vainu Bappu, Calar Alto, Teide, 
Foggy Bottom, Climcnhanga) participated only during the 
WEBT campaign in November 1997 (see Sect. 2.1) or 
around that period, while the others provided the long- 
term optical monitoring. All observers took frames by 
means of CCD cameras and performed data reduction us- 
ing standard procedures; the source magnitude was ob- 
tained by differential photometr y, ado pting the reference 
s tars c alibrated by Smith et al. ( 1985 ) and Fiorucci et al. 

(S)- 

Although the intense source variability makes a com- 
parison among the different datasets difficult to perform, 
a general agreement was found where the data overlap in 
time. 

Many large-amplitude outbursts are visible from Fig. 
1 ; in the B band the minimum and maximum magnitudes 
observed were 15.09 + 0.18 (JD = 2443907) and 20.5 + 0.05 
(JD = 2449692), respectively. In the R band the range of 
magnitudes spanned is from 14.03 ± 0.08 (JD = 2450811) 
to 19.47 ± 0.06 (JD = 2450278). 

An enlargement of the BVRI light curves in the period 
1996-2000 is shown in Fig. 2; an inspection of the i?-band 
curve, the best sampled one, reveals three major peaks, 
one for each of the last three observing seasons. 

Details of the 1997-1998 season are visible in Fig. 3, 
where different symbols refer to different observatories (see 
caption to the figure), with the exception of dark dots, 
which indicate data collected during the WEBT cam- 
paign. 



The peak at JD = 2450811 was observed by both the 
Perugia and the Roma Observatories, and corresponding 
peaks were seen also in the BVI bands (see Figs. 1 and 
2). What is quite impressive is the extreme sharpness of 
the flare: a brightness increase of 1 mag in 4 days (inside a 
2.8 mag rise in 27 days) was followed by a 1.3 mag decrease 
in only one day, leading to a 2.0 mag fall in 13 days. 

A smoother behaviour characterizes the 1998-1999 ob- 
servational season (see Fig. 4): the source was found to 
be very bright at the beginning of the season, skimming 
R = 14, and then a continuous decrease, interrupted by a 
couple of flares, led to R ~ 19, with a jump of 5 mag. 

In the last observing season the temporal coverage was 
definitely worse than in the previous two seasons (see Fig. 
2); however, a quite interesting feature is the fall of 1.6 mag 
in 48 hours occurred in the R band at the beginning of 
September 1999 (JD = 2451425.657-2451427.657). 

2.1. The November 1997 WEBT campaign 

The WEBT is an international collaboration among as- 
tronomers from all the world with the aim of organizing 
optical campaigns on blazars of specific interest (Mattox 
19994 Pggbt Villata et al. |2000[). These campaigns, last- 



ing from a few days to several weeks, are generally under- 
taken in concert with observations at other wavelengths, 
and can be triggered by the discovery of a flaring state of 
an object, usually in the optical band. 

This was indeed the case for the first-light WEBT 
campaign: the November 1997 observations of AO 
0235+16 were started after the detection of a con- 
siderable brightness increase (about 1 mag in a week) 
at the end of October 1997 (Webb et al. [1997D . An 
EGRET target of opportunity observation occurred be- 
tween November 3 and 11, 1997, but only an upper limit 
of 16 x 10 -8 photons cm -2 s _1 above 100 MeV could be in- 
ferred (Hartman, private communication). Observations 
by R XTE did not hnd a high X-ray flux (Webb et al. 
|2000D . 
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0235+164 (Aug 97 - Mar 98) 



M '4 f + * 

4 t + 




B750 8B00 
Julian Date - 2442000 



Fig. 3. Light curve of AO 0235+16 in the Cousins' R band 
from August 1997 to March 1998; data are from Vainu 
Bappu (asterisks), Abastumani (open circles), Roma (tri- 
angles), Perugia (diamonds), Torino (squares), and Teide 
(crosses); data from the WEBT campaign occurred from 
November 1 to 11, 1997 (JD = 2450754-2450764) are plot- 
ted as dark dots 



0235+164 (Jul 98 - Jan 99) 



9100 9150 
Julian Date - 2442000 



Fig. 4. Light curve of AO 0235+16 in the Cousins' R 
band from July 1998 to January 1999; data are from 
Abastumani (circles), Roma (triangles), Perugia (dia- 
monds), Torino (squares), Teide (crosses), and NOT (as- 
terisks) 



servational effort of many observatories around the world 
led to a satisfactory time coverage. 

In order to avoid spurious variations due to possible 
systematic effects among data from different observato- 
ries, in the following only flux changes observed by the 
same telescope will be considered. A brightness increase 
of 1.25 mag in two days was detected at the beginning 
of the campaign and noticeable flux oscillations were ob- 
served all the time. Variations of 0.2-0.4 mag in 5-7 hours 
were found in a number of cases; very impressive is the 
dimming of half a magnitude in about 5 hours detected 
in the night between November 3 and 4. Similar very fast 
variat ions have recently been observed by Romero et al. 
( |2000|) . 



2.2. Spectral behaviour in the optical band 

One item that is interesting to investigate when study- 
ing the optical variability of blazars is whether the flux 
variations are accompanied by spectral changes. In the 
new data presented in this paper we identified 47 optical 
spectra collecting data in at least three filters taken by 
the same telescope, with the requirement that the delay 
from the first frame to the last one is not greater than 
50 minutes. We discarded 7 spectra which contained large 
errors and analyzed the remaining 40. A selection of these 
"well behaved" spectra is presented in Fig. || in the log v 
- \og(yF„) plane, using different symbols for data taken 
at different telescopes. When taking into account that 
some discrepancies in the spectral shape can arise from not 
completely equal photometric systems used in the various 
observatories, no significant variations can be recognized 
when the flux changes. However, a deeper inspection re- 
veals that the spread corresponding to the B band is larger 
than that relative to the / one. This would suggest that 
the flux variations are of larger amplitude at the higher 
frequencies, as already observed in other blazars. 

All 40 "well behaved" optical spectra were then fitted 
by a classical power law: F u oc v a with a x 2 minimiz- 
ing procedure. The results are shown in Fig. ^, where the 
spectral index a is plotted as a function of the flux in the 
R band and a values with errors greater than 0.25 have 
been eliminated. One can notice that for high flux levels a 
presents an almost constant value of —2.8 — 2.7, while it 
tends to decrease with decreasing flux, although there is a 
not negligible spread. This means that the spectrum ba- 
sically steepens when the source gets fainter, a behaviour 
that is common to blazars. 



The light curve obtained in the first 11 days of 
November 1997 is plotted in Fig. 5. Data from each obser- 
vatory are indicated with a different symbol, according to 
Table 1 (Column 7). Weather was rather bad in those days 
in most of western Europe and North America, so that the 
temporal density of the curve is far from ideal, and far even 
from the density r eached in m ore recent WEBT campaigns 
(see Villata et al. 2000 , 2001 ); however, the common ob- 



3. Radio and optical fluxes: visual correlation 
analysis 

Radio light curves of AO 0235+16 are shown in Fig. || 
data at 22 and 37 GHz are from the Metsahovi Radio 
Observatory, those at 4.8, 8.0, and 14.5 GHz are from 
the University of Michigan Radio Astronomy Observatory 
(UMRAO). The observations at the Metsahovi Radio 
Observatory were made with the 13.7 m antenna using 
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-11.5 - 




14.70 14.80 
log v [Hz] 



Fig. 6. A selection of optical spectra from Roma (trian- 
gles), Torino (squares), and Teide (crosses); Julian Dates 
(—2442000) are indicated on the right; straight lines are 
drawn only to guide the eye through points of the same 
spectrum; no vertical shift is applied 



standard ON/ON techniques with DR 21 used as a cal- 
ibration source. Details on the observing procedure and 
data reduction can be found in Terasranta et al. (1992, 
19981). The data from UMRAO were taken with the 26 m 



paraboloid of the University of Michigan. A description of 
the observing and data reduction procedures is included 
in Aller et al. ( [1985| , |1999| ). 

As in the optical band, also in the radio AO 0235+16 
presents intense activity at all wavelengths, with pro- 
nounced outbursts lasting from several months to a few 
years. The overall flux variations (maximum value over 
minimum one) detected in the various bands are: 18, 22, 
32, 28, and 31 at 37, 22, 14.5, 8.0, and 4.8 GHz, respec- 
tively. In particular, in the last four years (JD > 2450300) 
variations up to a factor 18 were observed. 

In Fig. U radio fluxes (in Jy) are compared to the opti- 
cal ones (in mJy) , obtained in the following way: all mag- 
nitudes in the B band before JD = 2449000 have been 
transformed into R ones adopting the mean colour index 
< B — R > = 1.65 + 0.16 (which was derived by consider- 
ing all the B-R pairs from a same observatory separated 
by no more than half an hour); these data plus the real R 



F R [mJy] 

Fig. 7. The spectral index a as a function of the R flux 
for the 40 optical spectra from Roma (triangles), Torino 
(squares), and Teide (crosses) 



magnitudes after JD = 2449000 have been converted into 
fluxes by adopting Rieke & Lebovski (1985) and the law 
by Cardelli et al. ( |1989| ), and using a Galactic extinction 
A B = 0.341 mag (from NED). 

The existence of radio-optical correlations for AO 
0235+16 was investigated in a number of previous works; 
evidence for a simultaneous radio and optical variability 
was found in corresponde nce to the optical fl ares o f 1975 
and 1 979 (M acLeod et al. |l976|; Le dden et al. [l976| ; Rieke 
et al. |1976| ; Balonek & Dent |l980j ) and, more recently, to 
that which occurred in 1997 (Webb et al. |2000|). 



Clements et al. (1995) analyzed optical data taken in 
the period 1977-1991 and radio data at 8.0 GHz from 
UMRAO with the Discrete Correlation Function (DCF), 
and found that: "Overall, radio events lag optical events 
with lag times varying from to 2 months.". Takalo et al. 



(1998) visually compared optical data from 1980 to 1996 
with combined 22 and 37 GHz data from the Metsahovi 
Radio Observatory and noticed that some of the opti- 
cal spikes appear to be coincident with radio flares, while 
others have no counterparts. Moreover, the general trend 
looked very similar in both frequency regimes, suggesting 
some kind of correlation. 

A visual inspection of Fig. || shows that the big opti- 
cal outburst of 1975 has a big radio counterpart at 14.5 
and 8.0 GHz, while in 1979 a noticeable optical peak corre- 
sponds to a modest radio peak. Prominent radio outbursts 
at 22, 14.5, 8.0, and 4.8 GHz were observed in 1982, to- 
wards the end of the optical season, so that a possible 
optical peak might have been missed. A strong bright- 
ness increase was detected in 1987 in both optical and 
radio bands, but in this case the double-peaked optical 
flare seems to preceed the radio ones. The behaviour of 
the long radio outburst of 1990-1991 appears more com- 
plex; during the radio outburst a sharper optical flare was 
detected but not followed in details. The radio outburst in 
1992-1993 was double-peaked; just before the first radio 
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peak, an optical flare was detected; no other optical data 
were taken at the time when the radio fluxes reached their 
maxima. One interesting feature, however, is that if one 
looks at the better-sampled UMRAO data, the first radio 
peak seems delayed when proceeding from the higher to 
the lower radio frequencies. Indeed, the maximum value 
was reached on October 4, 1992 at 14.5 GHz, on October 
13 at 8.0 GHz, and on November 3 at 4.8 GHz. A radio 
flux increase in 1994 was practically not followed in the 
optical band. 

From these considerations it is clear that the main dif- 
ficulty in performing a meaningful study on possible radio- 
optical correlations is the paucity of optical data. 

The situation has been noticeably improved in the last 
years, because of the intense observational effort of the 
monitoring groups involved in the present work. Indeed, 
the big outbursts occurred at the end of 1997 and in 1998 
were accurately followed in all the radio bands and in the 
optical one (Fig. ^). 

A preliminary discussion on the radio-optical correla- 
tions in this period was presented by Villata et al. ( 1999| ). 

On December 28, 1997 (JD = 2450811.33), the i?-band 
flux reached 9.19 ± 0.68 mJy (R = 14.03); four days be- 
fore it was 3.90 ± 0.11 mJy and 24 hours after the peak 
the flux had dropped to 2.85 ± 0.13 mJy. This very sharp 
peak was observed by two different groups (Perugia and 
Roma) and also in different bands (see Figs. [I] and ^J). 
At the time of this optical flare, the radio fluxes were in 
a rising phase: they reached their maximum values later, 
likely first the shorter (observed by Metsahovi) and then 
the longer (UMRAO) wavelengths. The peaks at 37 and 
22 GHz were detected on March 3 and 2, 1998, respec- 
tively, although a gap in the radio datasets does not al- 
low to establish whether the peak of the outburst had al- 
ready occurred. The maximum value of the 14.5 GHz flux 
was seen on April 14, 1998 (JD = 2450918), and that at 
8.0 GHz on April 5, 1998, but in this case the peak proba- 
bly occurred later. This peak-delay effect has already been 
quoted above for the 1992-1993 outburst. The radio out- 
burst at 22 and 14.5 GHz is clearly double-peaked; this fea- 
ture is recognizable also in the 8.0 GHz data. It is notice- 
able that the second radio peak could reasonably be con- 
temporaneous at 22, 14.5, and 8.0 GHz, the best sampled 
bands, and corresponds to a second optical flare detected 
at the beginning of the observing season, in July-August 
1998. The solar conjunction period prevented to follow the 
rising phase of the optical outburst, so that the possibil- 
ity that the optical peak also in this case preceded the 
radio ones remains open. A second brighter optical peak 
detected on August 7, 1998 (JD = 2451032.63, R flux of 
8.94 ±0.16 mJy) occurred when the radio flux was in a de- 
creasing stage. Another important, sharp optical flare was 
finally detected on September 4, 1999 (JD = 2451425.66), 
which again has a radio counterpart, whose peak shows 
several weeks of delay. 

The above discussion demonstrates that, notwith- 
standing the great observational effort of the last years, 
we are still far from having the sufficient sampling to 



derive firm conclusions on the radio-optical correlations. 
We can only notice that, in general, when the observa- 
tional coverage is sufficiently good, a long time scale ra- 
dio flux increase corresponds to a short time scale optical 
brightness increase, whose peak may precede the radio 
one. Moreover, there are at least two cases (the 1992- 
1993 and 1998 double-peaked outbursts) where a progres- 
sive time delay in reaching the maximum value is observed 
when passing from the higher to the lower frequency radio 
fluxes. 

4. Statistical analysis 

In this section we apply the Discrete Correlation Function 
(DCF) analysis to the data shown in Fig. || in order to 
investigate the existence of characteristic time scales of 
variability and of optical-radio correlations. 

The DCF is a method specifically designed for un- 
evenly sampled datasets (Edelson & Krolik 1988; Hufnagcl 
& Bregman 1992| ) , which also allows an estimate of the ac- 
curacy of its results. 

Given two datasets and bj , one has first to combine 
all pairs, calculating the unbinned discrete correlations: 



UDCF ( 



(a,i - a)(bj - b) 



where a, b are the average values of the two datasets, and 
o a , at their standard deviations. The DCF is obtained by 
binning the UDCFy in time for each time lag r: 



DCF ( r ) = ^E UDCT V r ) 



M 

where M is the number of pairs cii, bj whose time lag A.^ = 
tj — ti is inside the r bin. Spurious correlations can be 
found, of the order of ±Af -1 / 2 . The standard error for 
each bin is: 



odcf(t) 



1 



M - 



UDCF, 



DCF(r)] 2 } 



1/2 



A positive peak of the DCF means correlation, which is 
stronger as the value of the peak approaches and exceeds 
one. A negative peak implies anticorrelation. Moreover, 
the width of the peak must be comparable to those of the 
autocorrelation functions, obtained by applying the DCF 
to each dataset coupled with itself. 

A preliminary binning of data in time before calcu- 
lating the DCF usually leads to better results, smoothing 
out flickering. The size of this binning is crucial especially 
in the optical, where short-term variations are frequent, 
since it can remove important information. Furthermore, 
an increase of the data binning interval implies an increase 
of the spurious correlations, while an increase of the DCF 
bin size has the opposite effect. Also the choice of the 
DCF binning is a delicate point, determining the balance 
between resolution and noise. In general, a similar value 
of M for each DCF bin and a limit of 10% to the appear- 
ance of spurious correlations must be assured in order to 
get reliable results. 



CM. Raiteri et al.: Optical and radio variability of AO 0235+16 



7 



Another method frequently used for searching charac- 
teristic time scales of variability is the Discrete Fourier 
Transform (DFT) spectral analysis for unevenly sampled 
data. We have adopted the implementation of the Lomb 
normalized periodogram method (Lomb 1976) discussed 
by Press et al. (1992). The presence of a sinusoidal com- 
ponent of frequency ujq = 2irv in the dataset is revealed 
by a large value of the periodogram P(lo) at to = ujq. 
The significance of the peaks is estimated by the false 
alarm probability, i.e. the probability that a peak is of 
height z or higher if the data are pure noise. It is given 
by F = 1 — (1 — c~ z ) M , where M is the number of in- 
dependent frequencies. Since we have scanned frequencies 
up to the Nyquist frequency v c = N / (2T) that the N data 
would have were they evenly spaced over the period T, we 
have set M — N. For very clumpy datasets as we have 
when considering the optical light curve, the value of M 
(and hence the false alarm probability) may be overesti- 
mated, and the significance of the peaks is consequently 
underestimated (see also Home & Baliunas 1986 ). 



4.1. Search for characteristic time scales of variability 

The optical DCF autocorrelation is plotted in Fig. [To] (top 
panel): as expected, it is dominated by flicker noise and 
reveals an important double-peaked signal at r ~ 4100- 
4200 days (11.2-11.5 years), and not negligible ones at 
r - 1200 and - 2100 days (~ 3.28 and - 5.75 years). 
These can thus be regarded as characteristic time scales 
of optical variability for AO 0235+16. A visual inspection 
of the optical light curve in Fig. || confirms these features. 

The peak at t ~ 5000 days in the optical autocorre- 
lation is just one point affected by a large error, deriving 
from the correlation between the 1979 and 1992 outbursts. 

The fact that the maxima at about ±4200 days exceed 
1 is due to the choice of a 2 day binning on the origi- 
nal dataset; a larger time interval for data binning (e.g. 
20 days) would reduce the importance of these features 
with respect to the central maximum, but would force the 
choice of a much larger DCF bin (~ 100 days) in order to 
avoid important spurious correlations. This in turn would 
imply missing important details. 

Figure [lO] also presents the DCF autocorrelation for 
UMRAO radio fluxes. The peaks are wider than in the 
optical case, reflecting the broader outbursts of the ra- 
dio light curves when compared with the optical ones. 
An important double-peaked maximum appears at 3900- 
4200 days (10.7-11.5 years) at all three frequencies, whose 
similarity with the optical one suggests radio-optical cor- 
relation. Other noticeable peaks are seen at 5300-5400 
and 6100-6200 days (14.5-14.8, 16.7-17.0 years) in the 
DCF autocorrelation function for the 8.0 GHz band. At 
14.5 GHz the former appears reduced, while the latter, 
deriving mainly from the coupling of the 1975 and 1992 
outbursts, is enhanced. At 4.8 GHz both disappear, but 
in this case the light curve is less extended in time; in 
particular, the 1975 outburst is totally missing. 



As for the Metsahovi radio light curves at 22 and 
37 GHz, their less dense sampling and more limited time 
extension lead to higher spurious effects, especially for 
large values of r. However, their autocorrelation functions 
confirm the signal centred at r ~ 4000 days, and present 
noticeable peaks at about 2000 days, that is the time sep- 
aration between the outbursts detected in 1987, 1992, and 
1998. 

The most interesting point emerging from the above 
discussion is the ~ 11.2 year characteristic time scale 
of variability, which is common to the optical and radio 
fluxes. By looking at the light curves in Fig. ||, one can 
understand this result by noticing that at a distance of 
about 4100 days the 1975 peak correlates with the 1987 
one, and this latter with the 1998 flare, while the 1982 out- 
burst correlates with the 1992 one. This means that there 
seems to be a ~ 4100 day characteristic variability time 
scale intersecting another ^4100 day time scale, which on 
one side is somehow surprising. On the other side, both 
the optical and the Metsahovi radio autocorrelations sug- 
gest that there may be a "periodicity" of about half the 
above time scale. Indeed, the 25 year time extension of 
the AO 0235+16 light curves would allow to interpret this 
halved time scale in terms of periodicity. The point is to 
understand why UMRAO radio autocorrelation functions 
do not show a strong signal at r ~ 2050 days. The reason 
is that this signal is damped by the delay of the 1982 out- 
burst. Indeed, by looking at the 8.0 GHz light curve (Fig. 
||), the best sampled one, one can recognize five large- 
amplitude outbursts peaking at JD - 2442000 = 723.60, 
3031.44, 4901.27, 6908.74, and 9000.00, spaced by 2308, 
1870, 2007, and 2091 days, respectively. The average pe- 
riod would thus be 2069+ 184 days, i.e. 5.67 + 0.50 years. 
The 8.0 GHz light curve folded assuming a period of 2069 
days is presented in Fig. O: the delay of the 1982 outburst 
is clearly visible. 

Notice that a weak signal at about 2000 days is ac- 
tually present in the 14.5 GHz autocorrelation function 
because at this frequency the 1982 outburst was preceded 
by a kind of pre-outburst that made the flux reach a high 
level earlier than the 8.0 GHz one. 

We have checked the reliability of the results obtained 
by the autocorrelation analysis by means of the Discrete 
Fourier Transform (DFT) technique for unevenly sampled 
data implemented by Press et al. (1992). In both the ra- 
dio and the optical cases we obtained many signals with 
significance levels F better than 0.001. In particular, the 
2069 day periodicity previously inferred is confirmed by 
the DFT analysis on UMRAO data. At 8.0 GHz, a clear 
maximum [P ~ 79) is found at frequencies v ~ 4.8- 
4.9 x 10 -4 day -1 , corresponding to periods of ~ 2050-2080 
days. Other strong signals in the 8.0 GHz data are found 
at 3.7, 2.8, and 1.8 years (P ~ 88, 46, 45, respectively). 
At 14.5 GHz the strongest maximum of the Lomb peri- 
odogram (P ~ 81) is right at v ~ 4.8-4.9 x 10 -4 day -1 , 
followed by the maxima (P ~ 68, 49, 48) corresponding to 
2.8, 1.8, and 3.7 year periods, confirming the results ob- 



8 



CM. Raiteri et al.: Optical and radio variability of AO 0235+16 



taincd for the 8.0 GHz data. Similar results are also found 
for the 4.8 GHz dataset. 

The DFT technique applied to the optical fluxes gives 
much more signals, making the spectral analysis rather 
complex. Surprisingly, the ~ 2100 day time scale discov- 
ered in the optical autocorrelation function and obtained 
by the DFT analysis of UMRAO data gives only a weak 
signal (significance level better than 3% only). To bet- 
ter visualize the matter, Fig. [l2| shows the optical light 
curve in the R band folded assuming a 2069 day period. 
Notice that such a period would not explain some impor- 
tant outbursts, in particular the major outburst observed 
in 1979. Other signals obtained by the DFT analysis at 
2.8 and 1.6-1.9 years confirm the time scales found for 
the radio fluxes, while there is not a strong optical signal 
corresponding to the 3.7 years in the radio. The strongest 
DFT signal corresponds to a ~ 200 day time scale, which 
is the time separation between the two peaks observed in 
the 1997-1998 outburst. 



4.2. Optical-radio cross-correlation 

The results of the DCF cross-correlation between data in 
the R band and the 8.0 GHz ones are shown in Fig. [l3] 
(bottom panel): the well-defined positive peak at r ~ 0-60 
days suggests optical-radio correlation, with optical varia- 
tions that can be both simultaneous and leading the radio 
ones by a couple of months. The DCF applied to the op- 
tical and 14.5 GHz datasets leads to a similar result (see 
Fig. [l^, top panel) : the radio variations appear correlated 
and delayed of 0-50 days with respect to the optical ones. 
These results are in agreement with what was derived by 
Clements et al. ( |1995| ). 

As for the radio-radio correlations, Clements et al. 
( |1995| ) found no time delay between the 14.5 and 8.0 GHz 
datasets and between the 8.0 and 4.8 GHz ones. However, 
as previously discussed, at least during the 1992-1993 and 
1998 outbursts, the radio light curves seem to indicate 
that the flux variation at the higher frequencies may have 
led that observed at the lower ones. Figure [l4| shows the 
results of the DCF cross-correlation between the 22 and 
8.0 GHz fluxes: in the top panel, where all data have been 
considered, the peak is not exactly centred at r = 0. This 
might suggest that the 22 GHz fluxes can lead the 8.0 GHz 
ones by several days. In the bottom panel of Fig. |l4| only 
data after JD = 2450300 were taken into account, so that 
only the 1998 outburst is considered: the delay effect in 
this case is enhanced. Thi s is in agreement with what was 
observed by O'Dell et al. (|l988l ) when analyzing the vari- 
ability of AO 0235+16 at eight radio frequencies, from 
318 MHz to 14.5 GHz. They found that flux-density vari- 
ations are clearly correlated, and events occur first at the 
higher frequencies and propagate to lower frequencies with 
decreasing amplitude. 
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Fig. 13. Discrete Correlation Function (DCF) between 
the R data and 14.5 GHz (top) or 8.0 GHz (bottom) data 
of AO 0235+16; all datasets have been binned over 2 days, 
while the DCF bin size is 20 days 



5. Discussion and conclusions 

Four years of optical and radio monitoring of the BL Lac 
object AO 0235+16 have confirmed the extreme variabil- 
ity for which this source is well known: a range of more 
than 5 mag was spanned in the optical band and a total 
variation of up to a factor 18 was registered in the radio 
fluxes. 

Optical spectral changes are not very evident, but our 
data suggest a spectrum steepening when the source gets 
fainter, a feature that has already been recognized for 
other blazars. This behaviour can be interpreted in terms 
of radiative losses in the electron population of the emit- 
ting jet: when the source is bright, cooling can be bal- 
anced or overcome by acceleration processes and the re- 
sulting spectrum is flatter; when the flux is low, radia- 
tive losses dominate and cause a spectral steepening, since 
the higher-energy ultra-relativistic electrons emitting syn- 
chrotron radiation cool faster than the lower-energy ones. 
However, also a geometrical interpr etatio n of the kind of 
that presen ted in Villata & Raiteri ( 1999 ; see also Villata 
et al. 200C) is possible: if a faint state of the source in a 
given band occurs when the portion of the curved jet emit- 
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Fig. 14. Discrete Correlation Function (DCF) between 
the 22 and 8.0 GHz data of AO 0235+16; in the bottom 
panel only data after JD = 2450300 have been consid- 
ered; both datasets have been binned over 2 days, while 
the DCF bin size is 20 days 



ting in that band becomes less aligned to the line of sight, 
then a spectral steepening is expected since the higher- 
frequency emitting portion of the jet departed first from 
the line of sight. 

Discrete autocorrelation function analysis of the opti- 
cal and radio light curves points out a characteristic time 
scale for the flux variability of AO 0235+16 of about 11.2 
years. However, a deeper insight into the optical and ra- 
dio light curves aided by DFT analysis and folded light 
curves reveals that major outbursts occur at roughly half 
the above time scale, even if one of the events (the 1982 
one) appears noticeably delayed. 

Indeed, in the best sampled 8.0 GHz light curve, there 
are five large-amplitude outbursts whose peaks are sepa- 
rated by 2308, 1870, 2007, and 2091 days, giving an aver- 
age "period" of 2069 + 184 days, i.e. 5.7 + 0.5 years. Taking 
into account that the AO 0235+16 redshift is z = 0.94, 
the relation P = P bs/(1 + z) implies that the period in 
the source rest frame is around 2.8 years. Other charac- 
teristic frequencies resulting from the DFT analysis of the 
radio light curves correspond to periods of 1.8, 2.8, and 
3.7 years. In the optical, the sparse sampling makes the 



analysis of variability time scales rather difficult. However, 
discrete autocorrelation function, DFT, and folded light 
curves are compatible with a ~ 2069 day periodicity. Such 
periodicity would not account for the observed big 1979 
optical outburst, whose radio counterpart was a rather 
modest flux increase. The optical autocorrelation analy- 
sis puts in evidence also another characteristic time scale 
of variability of about 1200 days (3.3 years), which corre- 
sponds to the time separation between the 1975 and 1979 
outbursts, while the DFT of the optical data confirm the 
1.8 and 2.8 year periods found in the radio data. 

A few previous studies investigated t he ex istence of 
periodicities in AO 0235+16: Webb et al. ( |l988[ ) analyzed 
the o ptical light curve with the Deeming DFT (Deeming 



1975 ), and derived periods o f 2.79 , 1.53, and 1.29 years. 



More recently, Webb et al. ( p000[ ) used unequal-interval 
Fourier transform and CLEAN techniques and obtained 
periods of 2.7 and 1.2 years for the optical variations. 
The time scales of long-term optical bas e-leve l fluctua- 
tions have been studied by Smith & Nair ( 1995 ) for three 
classes of AGNs; they found a best fit period of 2.9 years 
for AO 0235+16, attributing a moderate confidence to the 
estimate. The Jurkevic technique and the autocorrelation 
function were adopted by Fan (2001) for the analysis of the 
optical light curves, inferring periods of 1.56, 2.95 ± 0.25, 
and 5.87 ± 1.3 years. This last period is highly compati- 
ble with what is derived in the present paper. Roy et al. 
( |2000D made a cross-correlation analysis between optical 
and UMRAO radio data and applied Lomb-Scargle peri- 
odograms to the 14.5 and 8.0 GHz light curves to infer a 
~ 5.8 year periodicity. Their results are in agreement with 
what we find in the present paper. 

If the major outbursts in AO 0235+16 really occur 
every 5.7 ± 0.5 years, the next one should peak around 
February-March 2004, and a great observational effort 
should be undertaken since at least summer 2003 to get 
more information on the details of the flux behaviour 
in various bands. Unfortunately, the source is not visi- 
ble from ground-based optical observatories during spring- 
time because of solar conjunction. 

On the other hand, another great effort should be ad- 
dressed to the theoretical interpretation of such recurrent 
events, in particular to envisage possible mechanisms that 
can "disturb" the period. Many models have been pro- 
posed to explain the ~ 12 year period inferred from the 
optical light curve of another well-known BL Lac object, 
i.e. OJ 287; most of them foresee the existence of a bi- 



nary black hole system (Sillanpaa et al. 1988b; Lchto 
Valtonen |1996|; Villata et al. |1998j ; Valtaoja et al. |2000 



Abraham 2000] ). It would be interesting to see whether 



those models can be adapted to explain the quasi-periodic 
variability observed in AO 0235+16. 

As for the cross-correlation between the optical and ra- 
dio data, the analysis of the light curves is constrained by 
the sometimes very sparse optical sampling. However, in 
general we can say that a radio flux increase always corre- 
sponds to each major optical flare; the vice versa is likely 
true. Such correlation would suggest that the same mech- 
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anism is responsible for the optical and radio variations. 
However, two different behaviours seem to be present: op- 
tical and radio fluxes are seen sometimes to reach their 
maximum values at the same time, while in some cases the 
higher frequencies lead the lower ones. This would mean 
that two different mechanisms are presumably at work in 
the long-term variability. Outbursts occurring simultane- 
ously in all bands may suggest that the synchrotron emis- 
sion from the radio to the optical band is produced by the 
same population of relativistic electrons in the same re- 
gion. But they may also favour the microlensing scenario, 
which was already suggested for expl aining the vari ability 



of AO 0235 +16 ( Stickel et al. |1988j T akalo et al. |1998 ; 
Webb et al. 2000 ; but see Kayser 1988 for a critical dis- 
cussion). The presence of several foreground objects in the 
field of the source supports this scenario. One item against 
this hypothesis is that during the strong flux increase de- 
tected in both the optical and radio bands in 1997, the 
X-ray and 7-ray fluxes were not seen in a high state. 

On the other hand, the fact that in the radio and op- 
tical domains the variations observed at the higher fre- 
quencies lead those at the lower ones makes one think 
to an inhomogeneous jet, where a disturbance travelling 
downstream enhances first the emission at the higher fre- 
quencies and then the emission at the lower ones. An al- 
ternative view is given by geometrical models such as the 
previously mentioned helical model by Villata & Raiteri 
( |1999D , if we imagine that the portion of the jet emitting 
the higher-energy radiation gets closer to the line of sight 
before that producing the lower-energy flux, as already 
noted for the spectral changes. 

A different interpretation is required to explain the 
noticeable intraday variability shown by AO 0235+16 in 
both the optical and radio bands. Indeed, on short time 
scales, the flux behaviour can be very different with re- 
spect to that exhibited on long time scales. Kraus et al. 
( |l999| ) observed an "unusual" radio event in October 1992, 
in which the 20 cm maximum preceded the maxima at 3.6 
and 6 cm, and the variation amplitude was larger at the 
lower frequencies. They applied a number of models, both 
of intrinsic and of extrinsic nature, and concluded that, 
in any case, the size of the emitting region must be very 
small, implying a Doppler factor of order 100. Su ch a h igh 
value was subsequently confirmed by Frey et al. (200C) by 
analyzing VSOP observations of AO 0235+16. 

We found several episodes of very fast intraday vari- 
ability in the radio light curves; in some cases we are 
in the presence of single points that appear to deviate 
from the longer trends shown by the data, but there are 
cases in which the variation is confirmed by more than one 
point. Calculation of the brightness temperature for the 
fastest variability events (AF/Ai > 0.3 Jy day -1 ) in the 
14.5 GHz light curve (containing data generally of higher 
signal-to-noise ratio with respect to the lower frequencies 
ones) led to values by far exceeding the Compton limit, 
implying Doppler factors in the range 30-70. Indeed, such 
results must be taken with caution, since a much better 



sampling is needed in order to derive reliable Doppler fac- 
tors from the radio intraday variability. 

In conclusion, all the items pointed out by the present 
work (such as the 5.7 year periodicity, delayed flux vari- 
ations at lower frequencies, high Doppler factors) would 
need a further observational effort in the next years, in 
order to fix them to more precise results. 
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Fig. 1. Lon g-term monitori ng light curve s of AO 0235 +16 i n UBVRI; dat a are from Kinman & Rieke ( 1975 ), Rieke 
et al ( |1976|) , O'Dell et al (|l978a| , |l978b| ), P ica et al. ( |1980|) , Impey et al ( |l982| ), Barbi eri et al. ( |l982| ), Mole s et a l. 
( |1985| ), Smith et al. (jl987| ), Sillanpa a et al . ( |l988a| ), Webb & Smith ( |1989| ), Mead e t al. fll99C| ), Si tko fc Sitko ( [1991D , 
Sillan paa et al. ( |1991|) Takalo et al. fll992| ) Xie e t al. ( |1992[ ), Rabbette et al. ( |1996[ ), Webb et al. ( |1997| ), Takalo et al. 
( |'1998| ), Xie et al. ( 1999| ), and Ghosh et al. ( 200C ); data plotted after the vertical line are from the present work 
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Fig. 2. Light curves of AO 0235+16 in BVRI bands in the years 1996-2000; all data are from the present work 
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Fig. 5. Light curve of AO 0235+16 in the Cousins' R band during the first-light WEBT campaign of November 1997; 
for the explanation of symbols see Table 1 
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Fig. 8. Optical (mJy) and radio (Jy) light curves of AO 0235+16; radio data at 22 and 37 GHz are from the Metsahovi 
Radio Observatory, those at 14.5, 8.0, and 4.8 GHz are from UMRAO 
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Fig. 9. Optical (mJy) and radio (Jy) light curves of AO 0235+16 in the last years; radio data at 22 and 37 GHz are 
from the Metsahovi Radio Observatory, those at 14.5, 8.0, and 4.8 GHz are from UMRAO 
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Fig. 10. DCF autocorrelations for the optical and UMRAO radio fluxes shown in Fig. 
over 2 days, while the DCF was obtained with a bin size of 25 days 
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Fig. 11. The light curve of AO 0235+16 at 8.0 GHz folded assuming a period of 2069 days 
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Fig. 12. The light curve of AO 0235+16 in the R band folded assuming a period of 2069 days 



